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Abstract: We review recent research on the Milky Way galaxy and try to inves- 
tigate whether its shape is similar to other barred galaxies. The emphasis is given 
on microlensing research because this method can be useful in determining the 
shape of the Galaxy with the minimal set of assumptions. By analyzing plots of 
the microlensing optical depth, r as a function of galactic coordinates for different 
values of the axis ratio, q of the galactic halo, we have shown that observations are 
' best described by a flattened halo with 0.2< q <0.6. 

<n : 

■ Introduction 

Over the last decade a growing list of evidence led us to accept the hypoth- 
esis proposed back in 1964 by de Vaucouleurs (e.g. Kuijken 1996): our Galaxy is 
a barred one. The arguments were initially based upon morphological considera- 
tions. The development of observational techniques broadens the possibilities of 
investigation of the central parts of the Milky Way in several research fields. All of 
them, however, confirm de Vaucouleurs' idea that the Milky Way is a SAB(rs) type 
galaxy, having a bar, weak rings, and a four-arms spiral structure (Vallee 1995). 

We shall only mention here the following observational results (Kuijken 1996) 
that confirm these facts: 
, 1. photometric research that include: 

a) surface photometry (Blitz and Spergel 1991; Dwek et al. 1995); and 

b) star counts (Gould, Bahcall and Flynn 1997; Nikolaev and Weinberg 1997) 

2. kinematical research that include: 

a) gas kinematics (Binney et al. 1991); and 

b) stellar kinematics (Kuijken 1996, and references therein) 

3. gravitational microlensing (ML) researches. 

The existence of the bar successfully explains dynamical peculiarities observed 
in our Galaxy, and shows that they are normal features in a barred galaxy: 

- stellar and gas dynamics in the Galactic center region (Binney et al. 1991), 

- central activity definitely correlated with a presence of the bar in galaxies 
(Combes et al. 1995), 
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- a 2.6 x 10 6 M C3 black hole in the center of the Milky Way (Bower and Backer 
1998; Falcke et al. 1998), 

- the absence of HI and CO gas in the region 1.5 kpc < r < 3.5 kpc, implying 
that the major axis of the bar is r\, ar < r cr , where r cr = 2.4 kpc is the Milky Way's 
corotation radius (Binney and Tremaine 1987; Binney et al. 1991), 

- a molecular ring at 3.5 kpc, which can be explained by gas accumulation 
near the outer Lindblad resonance radius (^olr = 4.1 kpc in Milky Way) (Binney 
et al. 1991; Frcundrcich 1998), 

- asymmetry in the distribution of the red clump stars in the bulge (OGLE) 
(Stanek 1995), 

- asymmetries of the bulge photometric image (COBE-DIRBE) (Dwek 1995, 
Binney, Gerhard and Spergel 1996), and 

- excess of gravitational microlensing events compared to the theoretical esti- 
mates (Paczyhski et al. 1994) in the galactic bulge direction. 

All observations agree that the value of the bar inclination angle (to the Sun- 
Galactic Centre line) is between 10° and 30°. 

In this paper we will discuss gravitational ML research and try to establish the 
connection between the inner parts of our Galaxy (bar) and its outer parts (halo) 
using this new observational technique. The ultimate goal is to see how the Milky 
Way as barred galaxy can be compared to other barred spirals and whether some 
conclusions concerning the shape of the Galaxy can be drawn. 

Baryonic dark matter content of the Galaxy 

The dark matter (DM) content of the Milky Way is still unknown. From 
shape of its rotation curve (RC) (Mcrrifield 1992) one can see that a huge amount 
of mass still has to be identified. The difficulties in the determination of the RC 
led to uncertainties in the most important parameters such as the galactic constant 
Rq, which represents the distance to the Galactic center and the circular speed at 
the Solar radius, w (Merrifield 1992, Oiling and Merrificld 1998, Sackett 1997). 
Although the IAU 1986 standard values are: Ro = 8.5 kpc and v n — 220 kms -1 
some recent estimates allow the smaller values: Ro = 7.1 ± 0.4 kpc and vq — 
184 ± 8 kms -1 (Oiling and Merrifield 1998). In this paper we adopt the value 
Ro = 8. 5 ±0.5 kpc (Feast and Whitclock 1997) based upon an analysis of Hipparcos 
proper motion of 220 Galactic Cepheids and vo = 210 ±25 kms -1 that includes the 
best values from the HI analysis (v — 185 kms -1 ) and the estimated values based 
on the Sgr A* proper motion (v Q = 235 kms -1 ) (Sackett 1997). 

Without going into the discussions about the content of the DM in the halo, 
we only state here that one part (presumably smaller) has to be in the baryonic 
form. Namely, cosmic nucleosynthesis predicts that (Turner 1996): 

0.008 < n B h 2 < 0.024 (1) 

where CIb is the universal baryonic mass-density parameter (fi^ = 8irGps /3Hq) 
and 0.4< h <1.0. "Silent" h is used in parametrization of the Hubble constant 
i?o = lOO/i kms -1 Mpc -1 . Recent estimates (Fukugita, Hogan and Peebles 1998) 
give: 

0.007 < fl B < 0.041 (2) 
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with the "best guess" il B ~ 0.021 (for H = 70 kms" 1 Mpc" 1 ). 

Using the simplest dynamical estimate of the mass of the Galaxy (Kepler's 
third law): 

GM{r) = v{rf (3) 

where M(r) is the mass interior to r, v is the measured rotational velocity and r 
is the radius within which most of the light in galaxy is emitted. For luminous 
matter one can obtain: 

0.003< Q WM <0.007 (4) 

(Roulct and Mollcrach 1997, and references therein). This is consistent with severe 
limits on mass-to-light ratio in the Local Group imposed by deep blank sky surveys 
(Richstone et al. 1992; Hu et al. 1994; Flynn, Gould and Bahcall 1996), as well as 
with huge dynamical mass for the Milky Way inferred by Kulessa and Lynden-Bell 
(1992). 

The mass in the halo is dominated by the matter that is not, at least easily, 
detectable. So, one can write: 

^HALO>0.1>14r! LU M. (5) 

It can be seen the equations (2) and (4) that dark baryonic matter must exist; 
various types of such material have been suggested: gaseous clouds of plasma or 
neutral atoms and molecules, snowballs or icy bodies similar to comets, stars, 
planets, white dwarfs, neutron stars and stellar or primordial black holes (e.g. 
Peebles 1993). 

Microlensing methods and results 

In searches for the baryonic DM content the method of microlensing has so 
far proved successful. Its name derives from the fact that lensing of distant ob- 
jects is made by bodies with masses characteristic of a star or planet. Although 
the theoretical development of this idea started in 1964 (e.g. Peebles 1993, and 
references therein), it was the seminal paper by Paczyhski (1986) that showed that 
one can search for ML events in the Milky Way halo if it is made of stars or brown 
dwarfs. Rapid development of observational and computer technology led to the 
detection of a significant number of ML events (e.g. Mellicr, Bcrnardeau and Van 
Waerbeke 1998). The directions include Large and Small Magellanic Clouds (LMC 
and SMC) (Alcock et al. 1996, 1997b; Palanque-Delabrouille et al. 1998), Galactic 
bulge (Kiraga and Paczynski 1994) and M31 (Crotts 1996). 

All these surveys give results concerning two important parameters: masses 
of the intervening objects and the optical depth. In the Table 1 we give the tar- 
gets observed, names of the appropriate survey, mass ranges of the lenses, and 
corresponding optical depth. 
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Target 



Survey 



Mass range 



Optical depth 



LMC/SMC 


MACHO 


« 0.3 - 0.5 M Q 


tlmc 

TSMC 


= 2.9tJ-| x 10- 7 
= 1.5 - 3 x 10~ 7 


Gal. bulge 


MACHO;DUO;OGLE 


0.08 - 0.6 Mq 


Tbul = 


= 3.9±i"J x 10~ 6 


M31 


KPNO 


« 10 M Q 


TM31 


= 5 - 10 x 10~ 6 


LMC/SMC 


EROS2 


0.85 - 8.7 Mq 


TSMC 


= 3.3 x 10~ 7 



Table 1 . Targets in different ML surveys, the mass ranges of the lenses and optical 
depths. 

Another important quantity, the optical depth, r is used in discussion of ML 
and as we shall show, in determining the shape of the halo. It can be defined as 
the probability that at a given time a source star is being microlensed with an 
amplification larger than 1.34 (e.g. Roulet and Mollerach 1997). 

Here we wish to investigate in more details one property of the halo of the 
Milky Way that has often been neglected: its shape. It is known from the work of 
Sackett and Gould (1993) that instead of the equation for the mass density in a 
spherical halo: 



AttG 
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a 2 + r 2 



9(R T - r) 



(6) 



(where r is the Galactocentric radius, is the asymptotic circular speed of the 
halo, a is the core radius of the halo and Rt is the truncation radius) one should 
use the generalized formula: 



= tan^ z4 



V> 47rGVa 2 + C 



6(R T - () 



(7) 



where £ 2 



; tan 2 if) (z denotes height above the Galactic plane) . Here 



the flattening parameter ip is introduced: cos^ = q = c/a, i.e. its cosine is 
equal to the axis ratio and determines the shape of the halo En. En is related 
to q as q = 1 — n/10. Following Sackett and Gould (1993) we write the following 
expression for the estimate of the optical depth as a function of Galactic coordinates 
I (longitude) and b (latitude): 



r(l,b) = 



tan ip v 2 ^ 1 



dL(D - L)L 



tp c 2 D J Q (a 2 + R 2 ) - (2i?ocosZcos b)L + (1 + sin 2 b tan 2 4>)L 2 



(8) 

where we put Rq = 8.5 kpc and a = 5 (e.g. Alcock 1996). Now we integrate this 
equation and take D = 50 kpc (for LMC), D = 63 kpc (for SMC) and D = 770 
kpc for M31. Although Sackett and Gould (1993) take values for q starting with 
q = 0.4 (shape E6) we will start with admittedly extreme value q = 0.2 (shape E8) 
required by some theories such as DDM (decaying dark matter) theory (Sciama 
1997), based upon the recent Dehnen-Binney models of the Galaxy (Dehnen and 
Binney 1998). Attempts were made to show that this small value of q is not 
possible since q = 0.75 ± 0.25 (Oiling and Merrifield 1997), but at the cost that 
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Ro = 7.1 ± 0.4 kpc (Oiling and Merrifield 1998). We will nevertheless take into 
account such small value for q since we find DDM theory acceptable in solving 
different serious astrophysical and cosmological problems (e.g. Sciama 1993). 

There are several other lines of reasoning suggesting a high degree of halo 
flattening in spiral galaxies. One is for long time suspected (e.g. Ninkovic 1985) 
flattening of the Population II subsystem, which may be a consequence of the 
residual rotation, or more probably, global flattening of the gravitational potential 
created by dark matter. The other is the behavior of the gas distributed in the halo. 
If the seminal idea of Bahcall and Spitzer (1969) of extended gaseous halos of normal 
galaxies producing narrow absorption features in the spectra of background objects 
is correct, as indicated by recent low-redshift measurements (Bergeron and Boissc 
1991; Lanzetta et al. 1995), then the distribution of gas could tell us something 
about the shape of the gravitational potential. It is not a simple problem at 
all (see Barcons and Fabian 1987), but some results are quite suggestive. In an 
important recent paper, Rauch and Haehnelt (1995) have shown that for the most 
plausible values of Lya cloud parameters, the conclusion that their axial ratio 
(thickness/transverse length) is less than 0.25 is inescapable. This conclusion does 
not depend on the exact choice of model for Lya clouds, and, if the observations 
quoted above are correctly interpreted, would mean that the gaseous halos are 
also flattened by the same amount. One should keep in mind, though, that such 
absorption studies probe only "a tip of an iceberg" , since these objects are ionized 
to extremely high degree, and may as well contain dominant part of the baryonic 
density in eq. (1). 

Bearing this in mind, we solve the integral in the eq. (8) and give estimate for 
r in several cases of particular interest: 

• Optical depth r(Z, b) in the parametric space, with the parameter q fixed in 
steps of 0.2, i.e. q = 0.2, q = 0.4, q = 0.6, q = 0.8 and jwl. 

• Optical depth r for different targets: LMC, SMC, M31 and Galactic bulge 
(bar) in order to see what value of q determines the optical depth that is 
closest to observed value in the appropriate survey. 

Due to the space limitations, we hereby present just two three-dimensional 
plots. In the Figure 1, value of the optical depth r is plotted against galactic 
coordinates I and b. This is an estimate for q = 0.4, but can easily be done for 
other values. One can use such plots (on the same or smaller angular scales) in order 
to choose observing direction where the optical depth reaches maximal values. Such 
an example is shown in the Figure 2 where we plotted optical depth as function of 
coordinates I — 280. °5 and b = — 32. °9 of the Large Magellanic Cloud. Other plots 
and results of integration will be presented elsewhere. 1 

Results and conclusions 

After solving the integral in the eq. (8) for given values of the parameter q we 
looked for the values that match the optical depth obtained in various surveys. We 
found that the best agreement can be attained if we take 0.2< q <0.6. Namely: 



Some examples of plots and results of integration can be found in the postscript 



format at the following URL: http://www.geocities.com/CapeCanaveral/7102/ 



Belgrade-MACH0.html, or from the authors by e-mail. 
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1. For the case of the LMC, that has been studied rather well, the measured value 
of the optical depth based upon the sample of 8 events is r = 2.9^' \ x 10~ 7 
(Alcock 1997b) while we find that for q = 0.5 we have r«3x 1CP 7 (see Figure 
2). 

2. For the case of the SMC, that is studied less thoroughly, the optical depth is 
estimated as r = 1.5 — 3 x 10~ 7 (Alcock 1997c). Our results show that the 
model in the eq. (8) gives the value t>4 x 10~ 7 for q ps 0.5 and above. 

3. For the case of the galaxy M31 we found r ss 5 x 10~ 6 which is an accordance 
with the estimates 5 — 10 x 10~ 6 (Crotts 1996), under the assumption that 
9>0.2. 

4. Determining r towards the Galactic center is more complicated and we will 
not discuss it here. We only state that using the model in the eq. (8) we can 
estimate the halo contribution to the ML rate towards Galactic center which 
is between r « 5 x 10~ 8 (q — 0.6) and r w 1.6 X 10 -7 (q — 0.2); the estimated 
range for the total optical depth towards Galactic center is (r = 3.9^12 x 10 -6 ) 
(Alcock 1997a). 

From our estimates it is obvious that the spherical dark halo can be ruled out: 
the value for q lies in the interval: 0.2< q <0.6. Recent research shows that it is 
not uncommon case with spiral galaxies (Sackett and Sparke 1990; Sackett et al. 
1994; Oiling 1995). Very recently, the observations of the gravitational lens system 
B1600+434, consisting of two spiral galaxies (Gl and G2), where G2 is a barred 
one, suggest that it has g^0.4 (Koopmans, de Bruyn and Jackson 1998). 
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Fig. 1. Optical depth r as a function of the galactic coordinates I and b. Distance 
to sources is taken to be 50 kpc and the axes ratio is q = 0.4. 

Fig. 2. Optical depth r as a function of the galactic coordinates / and b for the 
Large Magellanic Cloud (LMC) for the axis ratio q = 0.5. 



8 



case q=0 . 4 




case q=0.5 LMC 



